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Abstract. This paper deals with the problem of hydrodynamic shear turbulence in non-magnetized Keplerian 
disks. Several papers have appeared recently on the subject, on possible linear instabilities which may be due to 
the presence of a stable stratification, or caused by deviations from cylindrical rotation. Here we wish to draw 
I attention to another route to hydrodynamic turbulence, which seems to be little known by the astrophysical 

. community, but which has been intensively discussed among fluid dynamicists during the past decade. In this 

' so-called bypass concept for the onset of turbulence, perturbations undergo a transient growth and if they have 

fT^ , initially finite amplitude they may reach an amplitude that is sufficiently large to allow positive feedback through 

■ nonlinear interactions. This transient growth is linear in nature, and thus it differs in principle from the well- 

known nonlinear instability. We describe the type of perturbations that according to this process are the most 
likely to lead to turbulence, namely non-axisymmetric vortex mode perturbations in the two dimensional limit. 
We show that the apparently inhibiting action of the Coriolis force on the dynamics of such vortical perturbations 
O [ is substantially diminished due to the pressure perturbations, contrary to current opinion. We stress the similarity 

of the turbulent processes in Keplerian disks and in Cartesian flows and conclude that the prevalent skepticism of 

^ ' the astrophysical community on the occurrence of hydrodynamic shear turbulence in such disks is not founded. 

■ 
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C3 1. Introduction by Balbus et al. (1996) and Balbus & Hawley (1998). 
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and who also insisted, as we shall do here, on the similarity 
of subcritical shear turbulence in plane and rotating flows. 

Our paper is organized as follows. In Sec. 2 wc out- 
line the fundamental problems of fluid mechanics bearing 
on the description of plane-parallel smooth shear flows. 
We present the recent developments in shear flow analy- 
sis made by the hydrodynamic community, and the novel 
bypass concept for the onset of turbulence in shear flows. 
In Sec. 3 we analyze the quadratic forms of the dynamical 
equations and discuss the essential role of pressure pertur- 
bations, which counteract the effect of the Coriolis force 
in 2D flows. In Sec. 4 we describe the transient growth 
of two dimensional perturbations. This will allow us to 
stress the similarity between the behavior of Keplerian and 
Cartesian flows, and the negligible role of the Coriolis force 
in the 2D perturbation kinematics and energetics. Finally 
we conclude in Sec. 5 that the bypass concept provides a 
plausible; sc;c'nario leading to turbulence in astrophysical 
disks. In the Appendix we present a simple sketch of the 
bypass scenario. 

2. From spectral decay to transient growth and 
bypass transition 

2.1. On the onset of turbulence in spectrally stable 
shear flows 

From the theoretical viewpoint there are flows that are 
spectrally stable at all Reynolds numbers (e.g. plane 
Couette or pipe Hagen-Poiscuillc flows), while some oth- 
ers become spectrally unstable at high enough Reynolds 
numbers (e.g. plane Poiseuille or Blasius flows). In the lat- 
ter case the flow is characterized by the critical Reynolds 
number that is the marginal value of this parameter above 
which the spectral instability occurs. In this sense the crit- 
ical Reynolds number of the spectrally stable flow is infin- 
ity. 

Here we shall consider smooth shear flows (without 
inflection point in the velocity profile), like in Keplerian 
disks. These are linearly stable according to classical fiuid 
mechanics. As shown by Rayleigh (1880), the existence 
of an inficction point (more precisely of a vorticity ex- 
tremum, as pointed out by Fj0rt0ft 1950) in the equi- 
librium velocity profile is a necessary condition for the 
occurrence of a linear (spectral) instability in hydrody- 
namic flows. Thus smooth shear flows (flows without a 
vorticity extremum) are "relaxed" in this context: they 
arc spectrally stable, meaning that exponentially growing 
solutions are absent. However, it is well known from lab- 
oratory experiments and from numerical simulations that 
finite amplitude perturbations may cause a transition from 
laminar to turbulent state at moderate, less then critical 
Reynolds number. 

This has led to the development of the concept of 
nonlinear instability in hydrodynamics (cf. Bayly 1986; 
Bayly et al. 1988; Herbert 1988; Orszag & Kells 1980; 
Orszag & Patera 1980, 1983). Until about ten years ago, 
the predominant view of this laminar-turbulent transi- 



tion was centered around the slow linear amplification 
of exponentially growing perturbations (the familiar T- 

5 waves), which modify the flow profile and thereby allow 
a secondary instability, further nonlinearity and finally a 
breakdown to turbulent fiow. According to this concept 
of nonlinear instability, the perturbations (and the turbu- 
lent state itself) are energetically sustained by nonlinear 
processes. This concept of nonlinear instability has been 
borrowed by astrophysicists, who still use it to explain 
turbulent processes in smooth astrophysical fiows, where 
no spectrally unstable solution is known, and in particu- 
lar in Keplerian disks fiows (see Balbus & Hawley 1998). 
However, there are subcritical transition phenomena that 
can not be attributed to the nonlinear instability. 

During the last decade of the 20th century, another 
viewpoint emerged in the hydrodynamic community on 
understanding the onset of turbulence in spectrally sta- 
ble shear flows, labeled as bypass transition (cf. Bobcrg & 
Brosa 1988; Butler & FarrcU 1992; Reddy & Henningson 
1993; Trefethen et al. 1993; Morkovin 1993; Gebhardt & 
Grossmann 1994; Henningson & Reddy 1994; Baggett et 
al. 1995; Waleffe 1997; Grossmann 2000; Reshotko 2001; 
Chagelishvih et al. 2002; Chapman 2002; Rempfer 2003). 
Although the bypass transition scenario involves nonlin- 
ear interactions - which intervene once the perturbations 
have reached flnite amplitude - the dominant mechanism 
leading to these large amplitudes appears to be linear. 
This concept of the onset of turbulence is based on the 
linear transient growth of vortex mode (aperiodic) per- 
turbations. The potential for the transient growth has 
been recognized for more then a century (see Kelvin 1887, 
Orr 1907a,b). However only recently has the importance 
of the phenomenon been better understood (cf. Moffatt 
1967; Marcus & Press 1977; Gustavsson & Huhgren 1980; 
Craik & Criminale 1986; Farrell & loannou 1993; Reddy 

6 Henningson 1993; Chagelishvih et al. 1997). 

The bypass concept implies that the perturbation en- 
ergy extracted from the basic flow by linear transient 
mechanisms causes the increase of the total perturbation 
energy during the transition process. The nonlinear terms 
are conservative and only redistribute the energy produced 
by the linear mechanisms. (A simple sketch of the bypass 
scenario in a Keplerian flow is given in the Appendix.) 

Thus, according to this concept, the transient growth 

of perturbations {i.e. the linear process) is the key element 
in the transition to turbulence in spectrally stable flows. 
The importance of this linear process has been stressed in 
the titles of several seminal papers: Henningson & Reddy 
(1994) ( "On the role of linear mechanisms in transition 
to turbulence"); Baggett et al. (1995) {"A mostly linear 
model of transition to turbulence") and Reshotko (2001) 
( "Transient growth: A factor in bypass transition"). 
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2.2. On the mechanism of perturbation amplification 
in spectrally stable shear flows 

The scheme of investigation of shear flow dynamics implies 
the following steps: introduction of perturbations into a 
mean flow, linearization of the governing equations, and 
description of the dynamics of perturbations and flow us- 
ing the solutions of the initial value problem {i.e. follow- 
ing temporal balances in the flow). In principle this can 
be done, but in practice it is a formidable task. Therefore 
the mathematical approach was changed. This was done 
by assuming that the solution is separable in eigenmodes, 
and then establishing the existence of at least one un- 
stable eigensolution. This approach became canonical in 
time, and resulted in a shift of ground: attention was di- 
rected to the asymptotic stability of the flow and no atten- 
tion was paid to any particular initial value or to the finite 
time period of the dynamics. Indeed, this phase of the evo- 
lution was not thought to have any significance - it was 
left to speculation. But recently the early transient period 
for the perturbations has been shown to reveal "rich" and 
complicated behavior leading to different consequences. 

It was found in 1990s that smooth shear flows are 
crowded by intense processes of mean flow energy extrac- 
tion by perturbations, energy exchange between pertur- 
bations, etc. even in the linear approximation, while fol- 
lowing the classical theory they are spectrally stable and 
consequently "relaxed." Especially it has been shown that 
a superposition of decaying normal modes may grow ini- 
tially, but will eventually decay as time goes on - a new 
linear transient channel of energy exchange between the 
mean flow and perturbations has appeared. Moreover, it 
has been shown that transient growth can be significant 
even for subcritical values of the Reynolds number and 
that its interplay with nonlinear processes can result in 
transition to turbulence without any "nonlinear instabil- 
ity" of the flows. 



2.3. On the nonorthogonality of linear operators and 
the nonmodal approach 

In fact an exact resonance is not necessary to obtain the 
transient growth of perturbation. This is the consequence 
of the non-normal character of operators that describe the 
linear dynamics of perturbations in flows (see, e. g. Reddy 
et al. 1993). The fact that the eigenfunctions of the lin- 
earized Navier-Stokes equations are not orthogonal {i.e. 
the operator is non-normal) is enough to allow for so- 
lutions that exhibit transient growth, depending on the 
initial conditions, before finally decaying (sec Criminalc 
& Drazin 1990; Trefethen et al. 1993). The mechanism of 
the transient growth is essentially inviscid - the operators 
are highly non-normal for large Reynolds numbers and the 
transient growth is asymptotically large in Re. (This fact 
is extremely important in the case of Keplerian accretion 
disks, where the Reynolds number is literally astronomi- 
cal: Re > 10i°!) 



These developments provoked a change of paradigm in 
the study of linear processes in the considered flows. It was 
the so-called nonmodal approach that became extensively 
used and even canonized in the 1990s on these grounds. As 
a result, substantial progress in the understanding of the 
shear flow phenomena has been achieved. The nonmodal 
analysis - some modiflcation of the initial value problem 
- implicates the change of independent variables from a 
laboratory to a moving frame and the study of temporal 
evolution of spatial Fourier harmonics (SFH) of perturba- 
tions without any spectral expansion in time (see Sec. 4). 

3. The similarity between Keplerian flow and 
plane parallel shear flow 

Since our aim is to elucidate the basic similarity of the 
dynamics of plane shear flows and of Keplerian disks, we 
approximate the Keplerian flow in the two dimensional 
(2D) limit by its tangent plane parallel flow, while re- 
taining the effect of rotation through the Coriolis force. 
This model, which ignores the purely geometrical compli- 
cations, is known as the shearing sheet model; it has been 
has been used for the study of the linear dynamics of both 
wave and vortex mode 2D perturbations (cf. Goldreich 
& Lyndcn-BcU 1965; Goldreich & Tremaine 1978; Drury 
1980; Nakagawa & Sekiya 1992 for wave mode pertur- 
bations, and Lominadze et al. 1988; Fridman 1989 and 
loannou & Kakouris 2001 for vortex mode perturbations). 

The dynamical equations are written in the local co- 
moving Cartesian co-ordinate system: 



(1) 



where (r, 0) are standard cylindrical co-ordinates and Q.q 
is the local angular velocity at r = tq: 



dVL X 
n{r) = 1^0 + -s-{r - ro) = f^o + A—. 

or ro 



(2) 



A = {rdr^)\r=ro ^^'■^ shcar parameter (in Keplerian 
disks A = -3rio/2 < 0). 

For the present purpose, we shall consider only 2D per- 
turbations, independent on the axial coordinate z. The 
resulting dynamical linear equations for the perturbations 
of the radial velocity (w^), the azimuthal velocity {u^) and 
the pressure (p) take the form: 



dUr dUr -„ dp 

-—- + Ax— 2iloU^ = -TT- + I^^Ur, 

at oy ox 

duA, ^ dus „^ . 9p ^ 

— — -I- Ax— h 2S2oWr -I- AUr = -TT + vI^ua,, 

ot ay ay 



dur du^ 
dx dy 



0. 



(3) 
(4) 
(5) 



The incompressible limit is taken here to leave out 
wave mode perturbations and to keep only the vortex 

mode perturbations, which are the basic ingredient of the 
bypass scenario. In this case the energy density of the 
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perturbation depends only on its kinematic characteris- 
tics and not on the thermodynamics, such as pressure. 

Pressure terms are also absent when the dynamical 
equations are written in vorticity form (for u = curlu), 
taking the curl of eqs. This is probably the reason 

why the pressure perturbations have been often ignored 
when discussing the dynamics of shear flows. However they 
play a very important role, even in the incompressible 
case, because they are the mediators of momentum ex- 
change between fluid particles, which results in the tran- 
sient growth of the vortex mode. (The physics of the "me- 
diator activity" of the pressure perturbations is described 
in detail in Chagelishvili et al. 1993, 1996.) The impor- 
tance of the pressure terms is clearly seen in the following 
analysis. 

We multiply eqs. and respectively by Ur and u^, 
in order to put them in kinetic energy form, and average 
them over a domain which is symmetrical in x. This pro- 
cedure is similar to that performed in papers by Balbus 
et al. (1996) and Balbus & Hawley (1998). 



dV 2 ' ^ ° 



■A] 



(6) 



(7) 



The source of instability is represented by the shear term 
A(urU(j,). But the Coriolis term 2riQ{urU^) is explicitly 
present too, and in eq. (|7|) it overbalances the source term, 
since A = — 3/2ilo in the Keplerian disk. It thus would 
seem at first sight that the Coriolis force has a profound 
influence on the shear flow stability, as it has been argued 
by Balbus & Hawley (1996). But this conclusion is contra- 
dicted by the fact that the Coriolis terms disappear when 
summing eqs. (j^ and |(7J), meaning that the growth rate 
of the total kinetic energy is independent of the Coriolis 
force. 

A more detailed analysis will show that the Coriolis 
terms may be eliminated also on the level of the dy- 
namical equations, by a suitable renormalization of the 
pressure perturbation. In this two-dimensional and incom- 
pressible case, the perturbation velocity field derives from 
the stream function ip: 

dip dip 



dy dx 
Renormalizing the pressure perturbation 



p 



eqs. OHHl El and [71 can be rewritten as follows: 



dur 
'dt 

du^ 
dt 

dUr 

dx 



A: 



dUr 

' dy 



dp" 
dx 



Ax—^ + AUr = 
dy 



dp" 
dy 



du^ 
dy 



0. 



(8) 
(9) 

(10) 

(11) 
(12) 



d ul 



dn" 

{Ur^)-iy{\VUr\'), 



(13) 



(14) 



In Sec. 2 we outlined the bypass concept in 
Cartesian/planar shear flows (a simple sketch of the by- 
pass scenario is presented in the Appendix) . Since we wish 
to apply a similar scenario to differentially rotating disks, 
we shall write the Cartesian counterparts of the averaged 
equations for comparison. If y is the streamwise variable, 
and x is the shearwise variable as in the considered case 
(i.e., X is analog to r, and y to cj)), then the equations 
are: 

dr 2' 



(15) 



dV 2' 



u{\VUy\^), 



(16) 



where A'^ is the shear parameter of the Cartesian flow 
Uq(Q,A'^x)] Ux, Uy and p'^ are the perturbations of ve- 
locity and pressure, respectively. If the flow is set up with 
U§ decreasing with x (i.e. A'^ < 0), then eq. lfTB|) shows 
that steady flow is marked by transport in the direction 
of increasing x, as in the Keplerian disk. 

These equations are identical to the normalized equa- 
tions derived above for the Keplerian case, which stresses 
the similarity of these two flows. But they deal with spa- 
tially averaged quantities, and one may wish to establish 
that the similarity is even more profound, namely that 
the local kinematics and energy of 2D vortex perturba- 
tions of plane and disk flows are the same in time and 
space, starting from identical initial conditions. 

4. Transient growth in 2D 

Mainly for the purpose of illustrating the bypass mech- 
anism, we shall describe now the evolution in time of an 
initial perturbation, in two dimensions. We perform a spa- 
tial Fourier transform of all relevant variables, as shown 
here for the pressure fluctuation: 



p(x, y, t) = p{kx{t),ky,t) exp [ikx{t)x + ikyy] 



(17) 



In a shearing flow, the perturbations cannot keep the form 
of a simple wave, since the wave-number of each spatial 
Fourier harmonics (SFH) depends on time (see Criminale 
& Drazin 1990 for a rigorous mathematical interpreta- 
tion). In our geometry 



k'xit) — kx(0) — Akyt, 



(18) 



thus the wave-number of each SFH varies in time along 
the flow shear: it "drifts" in k-space. Therefore p, in eq. 
(|17|l above, depends on time through kx{t). 

Equations (0), JHl and ((TJ take the following form in 
Fourier space: 



dt 



Ur{kx{t), ky, t)-2nu^{kx{t), ky,t) 
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= -ik,{t)p{k:,{t),ky,t), (19) 

^u<j>ik4t), ky, t) + {2n+A)urik^{t),ky, t) = 

^-ikyP{k,it),ky,t), (20) 

kx{t)ur{kx{t), ky,t) + kyU^{kx{t),ky, t) = 0. (21) 

For simplicity, we ignored here the viscous forces; their 
action may be easily incorporated in the analysis, as we 
shall see below. 

One readily shows that these equations possess a time 
invariant 



J = kyUr{t) - kx{t)u^{t), 



which expresses the conservation of vorticity in Fourier 
space. Making use of this invariant, we may write the so- 
lution of the system (|19II21|) as follows: 



Ur{kx{t^^ ky^ t) 



kqi 



u^{kx{t),ky,t) - -^jjjy^pj^^i 



p{kx{t), ky,t) 



i[Akl + 2no{kl{t)+kl)] 



X. 



(23) 
(24) 

(25) 



We see that the perturbed quantities vary aperiodically 
in time (as is natural of vortex mode perturbations) , and 
that they undergo transient amplification. Energy is ex- 
changed between the perturbations and the background 
flow, and that represents the basis of the bypass transi- 
tion scenario. It is important to note that the kinematic 
characteristics of the perturbations do not depend on the 
rotation rate JIq ~ they are identical to those of Cartesian 
shear flows. Only the pressure perturbation depends on 
JIq. Consequently, the energetics of the SFH is identical 
to that of the Cartesian case: 



u {kxit),ky,t) = u^{kx{t),ky,t)+u^{kx{t),ky,t) 

^U^{kx{0),ky,0)-'''- ' ^ 



kl(t) + kl 



kl{t) + kl 



(26) 



We see that the kinetic energy reaches maximum ampli- 
tude for kx{t) = 0. Similarly, the renormalized pressure 
perturbations do not involve rJo- 



P''{kx{t),ky,t)^p{kx{t),ky,t)-2^a^j{kx{t),ky,t) = 



lAkl 



[kl{t) + kiY 



(27) 



This expression is identical to that describing the evo- 
lution of the pressure perturbation in a Cartesian flow, 
which proves the similarity of transient growth in both 
types of flows, in two dimensions. 



The above analysis has been done in the inviscid 
case. Accounting for the viscous forces is a straightfor- 
ward procedure (see Fridman 1989): it may be done 
by multiplying the obtained equations by the factor 

V ^^dt'{kl{t') + ky) . For instance, eq. (O and 



exp 



will read as follows 

Ur{k'x ky . t) 



km 



kl 



' {kx {t)j ky , t) — 



km + ki 



-V /Jdt'[fc|(t')+fc^] 



2v llAt'{kl(t') + kl 



(28) 



(29) 



To estimate the maximum amplification which is 

(22) 

achieved in this transient growth at k^ ^ k\(^\ we as- 



sume that the largest perturbation wave-number corre- 
sponds to the Kolmogorov dissipation scale: |fca;(0)|ro « 
i?e^/^, where the Reynolds number is defined as Re = 
nor^/iy. Then 



kliO) 



Re 



3/2 



uHkx{0),ky,0) kl ' \{kyroy 
which can reach a huge value in astrophysical disks. 



(30) 



5. Discussion and conclusions 

We have seen that the concept of nonlinear instabil- 
ity in smooth spectrally stable shear flows has under- 
gone substantial revision in the hydrodynamic community. 
Nowadays the concept of the bypass transition to turbu- 
lence has become favorite and is under intensive develop- 
ment. 

Let us summarize the main features of this concept: 

- the onset of turbulence and the turbulent state itself in 
smooth spectrally stable shear flows is supported energeti- 
cally by the linear transient growth of vortex mode pertur- 
bations - the key ingredient of the turbulence is the vortex 
mode (eddy) perturbation, and the key phenomenon is the 
linear transient growth of the perturbation; 

- nonlinear processes do not contribute to any energy 
growth, but regenerate vortex mode perturbations that are 
able to extract shear flow energy. Doing so, nonlinear pro- 
cesses only indirectly favor the energy extraction by the 
vortex mode perturbations; 

- the non-orthogonal nature of the linearized Navier- 
Stokes equations is the formal basis of the transient 
growth; 

- the non-orthogonal nature increases with increasing 
Reynolds number; thus the operators are highly non- 
normal for the huge Reynolds numbers of Keplerian disks 
{Re > 10^'^) and the transient growth is asymptotically 
large in Re. 

In this paper we wanted to contribute to the revival 
of hydrodynamic shear turbulence as a possible expla- 
nation for the "anomalous" viscosity in non-magnetized 
Keplerian accretion disks. Specifically, we wished to draw 
the attention of astrophysicists to the bypass concept. 
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We showed that in disk flows there exist vortex mode 
perturbations that are similar to those which are held re- 
sponsible for the onset of turbulence in the Cartesian shear 
flow. The key point of our analysis is the interpretation of 
the important role of the pressure perturbations in the 
dynamical processes. 

In fact, the kinematics and energetics of the vortex 
mode perturbations are identical in the rotating disk and 
the Cartesian shear flows in the 2D case. The Coriolis force 
only causes deviation of the pressure perturbations from 
the Cartesian case. By focusing on the epicyclic motions, 
and underestimating the action of the pressure terms, 
one is led to the false conclusion that the Coriolis force 
suppresses hydrodynamic turbulence in Keplerian flows 
(Balbus et al. 1996; Balbus & Hawley 1998). 

This property of the pressure perturbations has been 
established here in the 2D case. In the more general case 
of 3D perturbations, the dynamics of vortex mode per- 
turbations is somewhat more complicated, as it will be 
shown in a forthcoming paper (Tevzadze et al.). But the 
role of the pressure perturbations is still to counteract the 
Coriolis force. In the 3D disk case, the wave-number do- 
main where the perturbation undergoes transient growth 
is smaller in comparison to 3D Cartesian flows, but this is 
compensated by the very large Reynolds number charac- 
terizing astrophysical disks. 

Other scenarios which may lead to hydrodynamic tur- 
bulence in astrophysical disks have been presented re- 
cently: they invoke a linear (spectral) instability arising 
from the stratification perpendicular to the disk (DubruUe 
et al. 2002) or due to deviations from cylindrical rotation 
(Urpin 2002; Klahr & Bodenheimer 2002). Such instabil- 
ities may well compete with the bypass mechanism pre- 
sented here, and it is not possible to conclude presently 
which is the best candidate for rendering astrophysical 
disks turbulent. 

In laboratory experiments the field narrows, because 
there is no stratification in the fiuid (other than imposed 
on purpose), and there the bypass mechanism provides an 
attractive explanation for the turbulence detected in flows 
which are linearly stable (angular momentum increasing 
outwards). We refer to Couette- Taylor experiments per- 
formed by Wendt (1933), Taylor (1936), Coles (1950) and 
Van Atta (1966). Very recently such turbulence has been 
observed also in rotation profiles which share the proper- 
ties of Keplerian disks, namely with their angular velocity 
decreasing outwards (Richard 2001). 

Decisive conclusions about the self-sustenance of the 
turbulent state need to be supported by numerical simu- 
lations. To our knowledge, the simulations reported so far 
have failed to detect hydrodynamic shear turbulence in ro- 
tating fiows with angular momentum increasing outwards. 
In a recent paper, Longaretti (2002) has discussed the pos- 
sible explanations of this negative result. The main reason 
is probably the lack of spatial resolution, which still pre- 
vents from reaching even the Reynolds numbers at which 
turbulence is detected in the laboratory. At this stage, we 
hope that the analysis presented here will help to build 



a fruitful background for the nonlinear numerical stabil- 
ity analysis of rotating astrophysical shear fiows. In any 
way, the definite answer will result from the convergence 
of theoretical, numerical and experimental work. 
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Appendix A: Bypass scenario of the transition to 
turbulence 

We present here a simple sketch of the bypass scenario 
applied to Keplerian disk flow, in the wave-number plane 
{kx, ky) (see Fig. lA.l| . We shall define as "active domain" 
the region where viscous dissipation may be neglected, i. e. 
where kl + k^ < kl, with k^, « i?e^/^/ro . 

The linear dynamics of the perturbation may be de- 
scribed by following each of its spatial Fourier harmonics 
(SFH) in this wave-number plane. We single out one that 
is located initially at some point 1 in the amplification 
part of the "active domain" , which meets the condition 
kx/ky < 0. According to eq. lfTH|) . as kx{t) varies in time, 
the SFH drifts in the direction marked by the arrows, since 
A < 0. (We present the drift of the SFH only in the up- 
per half-plane; since the perturbation is real, there is a 
counterpart in the lower half-plane.) Initially, as \kx(t)\ 
decreases, the energy of the SFH grows. This growth lasts 
until the wave vector crosses the line k^ = (point 2). 
Then, while undergoing attenuation, the SFH continues its 
drift until it reaches point 3, where it is dissipated through 




Fig. A.l. A simple sketch of the bypass scenario applied 
to Keplerian disk flow, in the wave-number plane (kx, ky). 
kx is the wave-number in the radial direction and ky in 
the azimuthal direction. The wave-number of the spatial 
Fourier harmonics (SFH) drifts from its initial position 
1, the SPH is amplified in quadrant II, reaches maximum 
amplitude in 2, is attenuated in quadrant I, and undergoes 
viscous dissipation in 3. But the amplification quadrant II 
is repopulated through nonlinear interaction from SFHs 
located in the attenuation quadrants I and HI. (See text 
for more details.) 

viscous friction. The same will occur with all other Fourier 
harmonics. Consequently, if the nonlinear interaction be- 
tween different Fourier harmonics is inefficient, the per- 
turbation disappears eventually. Permanent extraction of 
shear energy by the perturbations is necessary for their 
maintenance, which is possible when quadrants II and IV 
where kx/ky < are being repopulated through nonlin- 
ear interaction between Fourier harmonics of quadrants I 
and HI which have reached sufficient amplitude. This may 
be achieved through three-wave processes ki -t- k2 k; 
four wave processes ki -I- k2 -f ka ^ k; five wave pro- 
cesses, etc. In the figure we present an example of three 
wave process ki -I- k2 =J> k that contributes to the regen- 
eration of an SFH in the amplification area, transferring 
perturbation energy to it from the attenuation areas. The 
bypass scenario implies the dominance of this regeneration 
tendency of nonlinear processes, i.e. predominant transfer 
of perturbation energy from quadrants I and HI to quad- 
rants II and IV, in other words nonlinear positive feedback. 
To what extent the reproduction of the Fourier harmonics 
in quadrants II and IV is sustained, even in the case of 
positive feedback, depends both on the amplitude and on 
the spectrum of the initial perturbation. Nonlinear decay 
processes are weak at low amplitudes and are not able to 
compensate the linear drift of SFH in k-plane. As a re- 
sult, weak perturbations are damped without any trace, 
and without inducing transition to turbulence. The higher 
the amplitude of initial perturbation, and the stronger are 
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the nonlinear effects. At a certain amplitude (which, of 
course, depends on the initial vortex perturbation spec- 
trum in k-plane and on the Reynolds number), the non- 
linear processes are able to compensate the action of the 
linear drift, thus ensuring permanent return of SFH to the 
amplification areas (quadrants II and IV). This eventually 
ensures a permanent extraction of energy from the back- 
ground flow and maintenance of the perturbations, and 
thus of turbulence. 

Therefore the bypass scenario can be realized only in 
the case of finite amplitude perturbations and in each case 
it has a threshold that depends on the perturbation spec- 
trum and the Reynolds number. 

Let us summarize. According to the bypass concept, 
vortex mode (eddy) perturbations are the basic ingredient 
of hydrodynamic shear turbulence. The bypass scenario 
involves the interplay of four basic phenomena: 

- the linear "drift" of SFH in the k-plane; 

- the transient growth of SFH; 

- the usual viscous dissipation; 

nonlinear processes that close the feedback loop of the 
transition by mixing - by the angular redistribution of 
SFH in k-plane. 



